Flow instabilities in complex fluids: 
Nonlinear rheology and slow relaxations 
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Abstract. We here present two simplified models aimed at describing the long-term, irregular behaviours observed in the 
Theological response of certain complex fluids, such as periodic oscillations or chaotic-like variations. Both models exploit 
the idea of having a (non-linear) rheological equation, controlling the temporal evolution of the stress, where one of the 
participating variables (a "structural" variable) is subject to a distinct dynamics with a different relaxation time. The coupling 
between the two dynamics is a source of instability. 



Introduction 

Complex fluids are known to exhibit a wide range 
of unconventional behaviour when forced to flow, due 
to the intricate couplings between their structure at the 
mesoscopic scale and the imposed flow. In this work, we 
investigate theoretically certain situations where, under 
steady external drive, the long-term behaviour of the 
fluids is intrinsically unstable: the fluids never reach a 
steady state, and rather respond in a unsteady way (as 
shown by time measurements of e.g., the shear stress or 
the shear rate). Such situations include the appearance of 
sustained temporal oscillations in surfactant solutions 1 1 , 
01 as well as polymer solutions |3;]. In other cases, some 
apparently totally eiTatic temporal responses have been 
found, in dense colloidal suspensions |4] and several 
surfactant systems In the latter, there are some 

indications |5, 6] that the obtained signal is in fact the 
result of a. deterministic chaotic dynamics. Such a chaotic 
behaviour at virtually zero Reynolds number (and thus 
no inertia) must stem from a nonlinearity within the 
constitutive properties of the material, and has thus been 
dubbed "rheological chaos", or "rheochaos". 

General principles 

These unsteady behaviours raise numerous questions 
from the theoretical viewpoint. Various hypotheses can 
be made as to what mechanisms give rise to such tempo- 
ral instabilities (oscillations or rheochaos). The path that 
we are currently exploring relies on two main physical 
features shared by many complex fluids: 

(i) An underlying tendency to form shear-banded 
flows. 



(ii) A dependence of the present state of the fluid on 
past history, due to the presence of slow-relaxing 
structural modes. 

The appearance of shear-bands in shear flows is indeed 
one of the most frequent mechanical instability observed 
in complex fluids. Slow-relaxing "structural modes", on 
the other hand, are directly linked to the commonly ob- 
served memory effects - i.e., long-term metastability - in 
complex fluids, where the shear history of a sample can 
be remembered for a surprisingly long time. One way 
to interpret this is that the previous flow has disturbed 
in some way one or several quantities characterizing the 
structure of the fluid (for instance, the micellar length in 
a worm-like micelle system, or the local density in col- 
loidal systems, etc.), and that these structural variables 
then only relax over long periods, because they involve, 
e.g., collective motions, or are controlled by a slow, in- 
dependent physico-chemical process. 

In this article, we give a short overview of our ongo- 
ing work on this subject, presenting two related models 
which include the above-mentioned physical ingredients, 
as well as some preliminary results. 

A shear-thickening model with memory 

The first model is a toy rheological equation which 
describes a shear-thickening fluid with memory, within a 
purely scalar approach where only the shear component 
a of the stress tensor is considered. Imagining that the 
fluid is sheared in a Couette cell, we take z as the vorticity 
direction (the axial direction), and will consider spatial 
variations in that direction only. In the proposed model, 
the evolution with time t of the shear stress <7(z,f ) is then 



given by (in units where the elastic modulus is one) 

aiz,t) = 7-R{(T)-X j'M{t~t')a{z,t')dt' + KV^a{z,t) 

(1) 

with d = da/dt. The shear rate, 7, is supposed uniform 
in the z-direction (this is related to a low-Reynolds as- 
sumption). The positive term R{(y) ~ aa — bo^ + ca^ 
is a non-linear, instantaneous relaxation term, with a, b, 
c positive constants chosen so that ^(cr) has a decreas- 
ing portion. This creates a tendency for the fluid to form 
"vorticity" shear bands (stacked in the z-direction). The 
integral term over past states of the stress represents the 
delayed relaxation of slow structural modes, with A > a 
parameter (homogeneous to a jump rate) and M a decay- 
ing memory function; we will here specialise to an expo- 
nential memory, M{t) = T^"' exp(— j/Ts), and Ts will be 
the typical relaxation time of the slow structures. Finally, 
the inclusion of a non-local term, in the form of stress 
diffusion with diffusivity ff, is required to describe inter- 
faces between different bands in inhomogeneous flows. 

The original version of this model, proposed by Gates 
et al. |7], was purely temporal (no space variable z, 
JC = 0) and already proved capable of long-term unsteady 
responses: when 7 is externally fixed in a certain range, 
the stress displays sustained (space-homogeneous) oscil- 
lations akin to those of the van der Pol oscillator'. With 
the spatially-resolved version presented here, we are able 
to study the spatio-temporal dynamics of such unstable 
behaviour. 

Qualitative features of the model 

The qualitative features of the model are best ex- 
plained by rewriting the integro-differential equation {l) 
as an equivalent differential system^: 

a = y-R{G)-Xm+KV'^G , m^-- — - (2) 

The memory integral within eq. Q now appears as an 
auxiliary variable m{z,t), which follows its own dynam- 
ics and at each instant tries to relax towards the current 
value of the stress <7{z,t). 

The equilibrium flow curve (7(7) for the model is 
obtained when a = rii = 0. Flowing states on this curve 
are spatially uniform, and on each point of the curve, 
the memory m has had time to relax to the equilibrium 
value of the stress (m = a). However, when Ts is large 



Regimes of temporal chaos were also found when the relaxation 
term R was furthermore allowed to become slightly delayed. However, 
the physical meaning of this modification and its relevance remain 
uncertain. 

^ This transformation is possible with an exponential memory, but not 
in the general case. 
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FIGURE 1. Long-term flow curve (thick line) and underly- 
ing short-term curves (obtained, from left to right, with fixed 
values of the memory in = to m = 20, in steps of one). The 
stress range between the dotted lines is unstable. (Parameters: 
/f = 1, A = 40, Ts = 100, a = 100, b = 20, c = 1.02, K = 0.01.) 



(slow memory relaxation), this equilibrium flow curve 
may only be observed on long enough timescales. On 
timescales much shorter than Ts, the fluid will behave 
as if the memory m were static: one has therefore also 
to consider a set of "instantaneous" flow curves which 
correspond to the relaxation of the stress (i.e., & = 0) at 
fixed m. As the memory slowly evolves, the fluid will 
then accordingly "jump" from one instantaneous flow 
curve to the other. 

It is then easy to understand qualitatively why the fluid 
has an unstable behaviour in a certain range of stress. 
In Figure 1, the long-term, equilibrium flow curve is 
drawn together with the set of instantaneous curves: the 
sole inspection of the equilibrium curve, which has no 
decreasing portion (for the choice of parameter values 
considered), could let one think that all states are stable; 
but one observes that there is a region where the long- 
term curve is in fact crossed by decreasing parts of short- 
term curves. Consequently, at these intersection points, 
the fluid has an instantaneous tendency to destabilize, 
thereby precluding the establishment of the equilibrium 
state 12]. 

Numerical study of the model 

We now briefly present some of the results that have 
been found so far on the spatio-temporal dynamics of the 
model. 

The model was studied numerically through a spec- 
tral Galerkin truncation fs*], where cj(z,f ) and m{z,t) are 
decomposed as a finite sum of spatial Fourier modes, of 
the form (J„(f)cos(q'„z) with wave-vector ^„ —nn/H for 
the stress (H is the total height of our imaginary Couette 
device), plus a similar set of memory modes m„{t). (We 
usually take the first ten modes into account in our nu- 
merics.) The system of equations then reduces to a 



set of coupled ordinary differential equations governing 
the temporal evolution of (7„(f) and mn{t). 

As is the case for conventional shear bands, two differ- 
ent protocols can be followed in the numerical study of 
the instability (details will appear elsewhere |91): either 
working at fixed shear rate 7 (i.e., with a shear-controlled 
Couette cell), or working at fixed torque, or equivalently 
at fixed average stress (a) (stress-controlled Couette). 
These two protocols lead to rather different results, as 
we shall now see. 



Working at fixed y 

We here work with a fixed, externally imposed value 
of 7, and compute the values of the different modes 
cr„(f) and m„{t), from which we reconstruct C7(z,f) and 
m{z,t). For numerical solutions carried within the unsta- 
ble window, our results show that, in the vast majority 
of cases, the spatio-temporal dynamics of the instability 
is in fact essentially purely temporal: regardless of their 
initial magnitude, all the modes rapidly decay and disap- 
pear, except for the spatially uniform mode CJo(f)- This 
mode then oscillates alone, with a fixed period and a reg- 
ular shape, as in the temporal version of ref. fl^. 

In a few, small regions of the parameter space ex- 
plored, it has been possible to obtain a slightly richer dy- 
namics, where some of the lower modes survive and un- 
dergo small-scale periodic oscillations, alongside a large 
uniform mode oscillation which still dominates. 

Working at fixed (cr) 

We now work at imposed average stress (a); in terms 
of Fourier modes, this corresponds to fixing the value of 
the uniform Fourier mode ao{t) — const., which will thus 
not be able to oscillate at all. A very rich spatio-temporal 
behaviour then arises, as the higher spatial modes will 
now, by necessity, be involved in the instability. 

Different unstable features may appear In Fig- 
ure 2-a, we present a most striking one, called "flip-flop 
shear-bands": a shear-banded profile appears, with a low 
shear band and a high shear band separated by an inter- 
face; but these bands are unstable, and periodically flip, 
the higher band becoming the lower and vice-versa. The 
"flipping time" is very short as compared to the latency 
time between two flips. Figure 2-b shows the correspond- 
ing temporal evolution of the total stress C7(z,f ). 

Much more irregular-looking time variations can also 
be obtained for different parameter choices, as Fig. 2-c 
shows for the stress. We emphasize however that these 
patterns remain time-periodic; so far no chaotic response 
has been found in the present model. 
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FIGURE 2. (a) "Flip-flop" shear-bands for (cr) = 7.0: From 
top to bottom, successive snapshots of the stress (j{z,t) vs. 
position z, at times t = 8800,8840,8880,8920,8960. {b) Time 
series of the stress at position z = for the same value of (cr). 
{c) Time series of the stress at position z = for (a) = 9.0. 
The time window corresponds to two periods of the signal. 
(Parameters same as Fig. 1.) 

A model with "fluidity" 

We would now like to introduce another fluid model, 
based on the same general ideas as previously, but more 
specifically oriented towards solutions of wormlike mi- 
celles and polymers. The shear stress E in the fluid is 
locally the sum of a polymer or micellar part, a, and a 
Newtonian part tj 7 corresponding to the solvent: 

r = (7 + T77 (3) 

The model's equations are as follows: 

a = + —^g{YT) + KV^-a (4) 

T T 

The stress evolution equation @ has a classical form, 
and describes the relaxation (with a timescale t) of the 
"polymer" stress towards a value which, at a given 7, 
is controlled by G(T)g(7T). G(t) is the elastic modulus 
and generally depends on the value of T (see below). The 
function g is hump-shaped, hence conferring on the fluid 
a tendency to shear-thinning, and shear-banding, with 
bands in the velocity gradient direction (i.e., the radial 
direction of the Couette cell). 

Similarly to the model with memory, one of the vari- 
ables involved in the stress equation is subject to a dis- 
tinct dynamics: this is now the Maxwell time T of the 
fluid, which, as stated by eq. Q, relaxes towards a shear- 
rate dependent equilibrium value Teq(7) with a character- 
istic time Ts . We note that having a dynamical Maxwell 



time is indeed very similar to the "fluidity" model in- 
troduced by Derec et al. in the context of paste flow 
(see fld\). 

Here again, T is a "structural" variable in the sense that 
it reflect variations in the local structure of the fluid: for 
semi-dilute solutions of wormlike micelles, it will typi- 
cally relate to local variations of the mean chain length of 
the micelles. In this interpretation, the chain length dis- 
tribution, and consequently the mean chain length, may 
change only through the action of some micelle-micelle 
chemical reactions, which will have their own timescale 
Ts - this time will thus control the relaxation of T, as in 
eq. (|5}- On the other hand, the dependence of the equi- 
librium value Teq on the shear rate can be interpreted as 
a displacement of the chemical equilibrium by the me- 
chanical shearing (for example, decreasing Teq, by help- 
ing chain scission, or, increasing Teq, by helping poly- 
merisation). For more dilute solutions (non-entangled), 
one may rather interpret the variations of Teq as related to 
the strong shear-thickening transition usually observed 
upon varying 7: the structural changes occurring at the 
transition (gelation) will then strongly alter Teq. 

In accord with the model with memory of the previous 
section, we have focussed on slow relaxations of the 
structural variable: Ts ^ 1 . Then one can again construct 
a "long-term", master flow curve where all variables are 
equilibrated (d — t — 0), and an underlying set of "short- 
term" curves where t is fixed (quasi-static) and the sti'ess 
is equilibrated with respect to that value of T. As seen in 
the previous model, an increasing portion of the "long- 
term" curve is here also made unstable when crossed by 
decreasing portions of underlying curves. 

Two variants of the model can be studied, depending 
on whether the elastic modulus G is affected or not by 
changes in the structure of the fluid, that is, changes in 
T. Variant 1 corresponds to G(t) = Go being constant, 
which would be suitable for a semi-dilute, entangled 
solution of micelles, where (reasonable) changes in the 
mean micellar length leave the modulus unchanged. In 
variant 2, G varies with T, as would e.g. be the case in the 
shear-thickening transition of dilute micelles, where the 
significant structural changes affecting T will affect G as 
well. Work on both variants of the model is currently in 
progress, but results so far are promising. Figure 3 shows 
an example of results in variant 2. 

To conclude, we would like to mention extremely 
interesting results by S. Fielding and P. Olmsted fTTIl : 
working independently on essentially the same model 
as eqs. 0-(|5}, they have been able, within variant 1, 
to obtain chaotic-like signals in a regime where both 
the underlying curves and the master flow curve have 
decreasing portions, and where the structural relaxation 
occurs on timescales comparable to the stress relaxation 

(Ts 0^ 1). 

Many questions are still open in the study of the type 
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FIGURE 3. Upper figure: Long-term flow curve (thick line) 
and underlying short-term curves in variant 2. Lower figure: 
very irregular (though periodic) time series of the total stress 
for (y) = 4.0. 

of models which have been described in this article. Are 
the existence of distinct "structural" relaxations really 
the origin of the instabilities observed in experiments? 
Is the regime just described the only one displaying 
chaotic-like behaviour, and relatedly, what is required for 
regular or irregular-looking periodic motion to destabi- 
lize into chaos? 
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